Generating a Fractal Butterfly Floquet Spectrum in a Class of Driven SU(2) Systems 
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A scheme for generating a fractal butterfly Floquet spectrum, first proposed by Wang and Gong 
[Phys. Rev. A 77, 031405(R) (2008)], is extended to driven SU(2) systems such as a driven two- 
mode Bose-Einstein condensate. A new class of driven systems without a link with the Harper 
model context is shown to have an intriguing butterfly Floquet spectrum. The found butterfly 
spectrum shows remarkable deviations from the known Hosftadter's butterfly. In addition, the level 
crossings between Floquet states of the same parity and between Floquet states of different parities 
are studied and highlighted. The results are relevant to studies of fractal statistics, quantum chaos, 
coherent destruction of tunneling, as well as the validity of mean-fleld descriptions of Bose-Einstein 
condensates. 
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INTRODUCTION 



Hofstadter's butterfly spectrum of the Harper model 
first discovered in two-dimensional electron systems 
subject to a square lattice potential and a perpendicular 
magnetic field, has attracted tremendous mathematical, 
theoretical and experimental interests. For an arbitrary 
irrational value of one system parameter, the spectrum 
of the Harper model is a fractal, which has been strictly 
proved after decades of research on the "Ten Martini 
problem" Q]- As one important property of Hofstadter's 
butterfiy spectrum, the number of its sub-bands depends 
on the arithmetic property of the fiux of the magnetic 
field. As the magnetic flux changes smoothly and thus 
varies between irrational or rational numbers, the gap 
between the sub-bands shows fractal properties and will 
close itself infinitely times [3|. This implies that Hof- 
stadter's butterfly spectrum contains infinite quantum 
phase transition points 4]. 

Early quantum chaos studies established that the Flo- 
quet (quasi-energy) spectrum of periodically driven sys- 
tems may display a fractal butterfly pattern as well 
[1)[1)0|- However, the nature of the fractal Floquet spec- 
trum is still poorly understood for three reasons. First, 
because the eigen-phase of Floquet states is restricted to 
a range of 2n, understanding a Floquet spectrum asso- 
ciated with an inflnite-dimensional Hilbert space is sub- 
tle and challenging [8]. Second, a rigorous mathematical 
proof about the fractal nature of a butterfly Floquet spec- 
trum is still lacking. Third, previous flndings regarding 
to fractal Floquet spectrum were largely limited to the 
so-called kicked-Harper model (a driven version of the 
Harper model) [1, M, El and its variant [1, 0, [HI, [11 ■ 



'Electronic address: |phygj@nus.edu.sg| 



Because of great interests in studies of quantum con- 
tro l, espec ially in studies of dressed matter waves 

there are now promising possibilities for the 
engineering and simulation of driven ultracold systems 
with a prescribed Floquet spectrum. Reference [6| rep- 
resents a recent attempt in this direction. In particular, 
in Ref. we showed that by designing two (5-kicking 
sequences the Floquet spectrum of a double-kicked-rotor 
system can be made to be a Hofstadter's butterfly, and 
the spectrum is identical with that of a kicked-Harper 
model when a certain parameter takes an arbitrary irra- 
tional value (l3 |. 

In this paper, we reveal yet another class of butter- 
fly Floquet spectrum using driven SU(2) systems, which 
are realizable by, for example, a driven two-mode Bose- 
Einstein condensate (BEG). As seen below, the basic 
strategy is essentially an extension of our previous work 
[1] , thus suggesting the possibility of synthesizing butter- 
fly spectrum in many other systems that go well beyond 
the context of two-dimensional electron systems or the 
Harper- model context. 

Our findings about the butterfiy spectrum of driven 
SU(2) systems are both motivating and fascinating. First 
of all, as explained below, now all the three popular 
paradigms of quantum chaos, i.e. the kicked-rotor model, 
the kicked-Harper model, and the kicked-top model, are 
linked together, insofar as any one of them can be used to 
generate quantum critical systems with fractal statistics. 
Second, the butterfiy spectrum obtained in driven SU(2) 
systems is significantly different from Hofstadter's but- 
terfiy, with remarkable aspects. For example, we show 
that with one certain system parameter fixed the overall 
butterfly pattern is insensitive to the number of bosons 
(denoted N) in the BEG, but some detailed features de- 
pend on whether N is odd or even. This may serve as 
a useful guide for seeking dramatic coherence effects in 
a BEG. Another interesting aspect is that the butterffy 
spectrum contains many level crossings between states of 
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different parities and thus many points of coherent de- 
struction of tunneling (CDT) [l8[, with the total number 
of CDT points found to scale as ~ iV'^". As an analog of 
quantum phase transitions in driven systems, the found 
butterfly pattern also contains many level crossings be- 
tween same-parity eigenstates. Due to these distinctive 
properties, the butterfly spectrum reported here repre- 
sents a new class of fractal Floquet spectrum and may 
become a test bed for a number of research topics. To em- 
phasize the remarkable differences between the butterfly 
spectrum found here and that associated with Harper's 
model, we refer to the newly found spectrum as "butterfly 
spectrum" instead of "Hofstadter's butterfly spectrum" . 

The main results of this study has been briefly reported 
in Ref. [l^ and this paper represents a full-length de- 
scription of our flndings. In Sec. II we will introduce the 
model SU(2) system and explain the main idea behind 
our study. In Sec. Ill we study the peculiar multifrac- 
tal spectral properties of the butterfly spectrum and the 
associated level crossings. The relevance of the underly- 
ing classical limit is also discussed in detail. To motivate 
possible experiments, we discuss some related issues in 
Sec. IV. We conclude this study in Sec. V. Appendices 
A, B, and C present some further details that may be of 
interest to some readers. 



II. DRIVEN SU(2) MODEL 

Our driven SU(2) model was motivated by a driven 
two-mode BEC system, proposed earlie r fisl [20l [2l| to 
realize the well-known kicked-top model [24I in the quan- 
tum chaos literature. In a very general form, a driven 
two-mode Bose-Hubbard model can be written as 

H = f{t)h{a\a2 + ajfli) -I- g{t)Ti{a\ai - 0^02)^, (1) 

where and a\ are the bosonic annihilation and cre- 
ation operators for the ith mode, f{t) describes the time- 
dependent tunnelling rate between the two modes, and 
the g{t) term describes the self-interaction between same- 
site bosons, whose time dependence can be achieved by 
Feshbach resonance induced by an additional magnetic 
field. Note that the total number of bosons N = a\ai + 
a\a2 is a conserved quantity. For a fixed N , the dimen- 
sion of the Hilbert space is + 1. Using the Schwinger 
representation of angular momentum operators, namely, 
Jx — (0102 + a2'^i)/2, Jy — {a\ai — a|a2)/(2i), and 
Jz — {a\ai — 0202)72, Eq. ^ reduces to 

H = 2f{t)hJ,+Ag{t)hJl (2) 

This above Hamiltonian makes it clear that its dynam- 
ics is solely determined by the SU(2) generators Jx, Jy 
and Jz. The total angular- momentum quantum num- 
ber J is given by J = N/2. The Hilbert space can be 
expanded by the eigenstates of Jz, denoted |m), with 
Jz\fn) — m\m). The population difference between the 



two modes is given by the expectation value of 2Jz. It is 
also important to note that if we exchange the indices of 
the two modes, then J^ is invariant, Jz — *■ —Jz, and as 
a result the Hamiltonian in Eq. ([2]) is unchanged. This 
reflects a parity symmetry of our model, which will be 
exploited below. 

Consider then two specific forms of f{t) and g{t). In 
the first case f{t) = a/(2T), g{t) = goJ2n{^i^ ~ ■^''^^ ~ 
t) — 6{t — 2nT)]. The Floquet operator, i.e., the unitary 
evolution operator F from 2nT + 0+ to {2n + 2)r + 0+, 
is then given by 

F — e'''''^^/(^-^)e~*"''="e~*''''^^'^'''e^'"'^^, (3) 

where rj — Ag^N . Interestingly, the first two or the 
last two factors in Eq. ([3]) constitute the Floquet op- 
erator for a standard kicked-top model [22]. As such 
our driven system here can be regarded as a "double- 
kicked-top model". Alternatively, if we set g{t) = go/^,, 
f{t) — ^ J2ni^i^~''^^)~^^i^~^'''~0]^ where ^ is the time 
delay between the two delta kicking sequences, then the 
associated propagator F' from nr — 0+ to {n + 1)t — 0+ 
is given by 

Under the special condition 4gor/^ = 2fc7r (Skir) for inte- 
ger J (half integer J), where k is an integer , the factor 
g-i(4goT/«)Jf yj^i^y ^jjg (2J+ l)-dimensional Hilbert 
space and hence F' becomes identical with F. Based on 
this, one now has two different scenarios for realizing F, 
the key operator to be analyzed below. 

To explain our motivation of considering the F op- 
erator, let us consider the \m) representation. In that 
representation the third factor e^"''^^/'^'^'' of F equals 
g-ji7m /(2j)^ which is a pseudo-random number for irra- 
tional rj/J. Interestingly, the first factor of F however 
effectively induces a time-reversal of the third factor and 
thus partially cancels this pseudo-random phase. Indeed, 
using the SU(2) algebra |23|, the product of the first three 
factors of F in Eq. ^ is given by 

g»»7jf/(2J)g-iaJ,g-ir,jf/(2J) 
^ g-«a{(./x/2+iJ„/2)e-[''<2'^-+i)/<2'^)l+c.c.} j-g-j 

This shows that the 77-dependent term entering into F be- 
comes e*!'''^'^^"'"^-'/^^'^''!, which is always a quasi-periodic 
number e^^''^^'""''^^/''^'''-'! in the |m) representation. Ac- 
cording to our early work [11], such a partial cancela- 
tion of quasi-random dynamical phases implies intriguing 
spectral properties. 

For later discussion we also study the classical limit of 
F. To that end we consider scaled variables x = Jx/J, 
y — Jy/J^ z — Jz/J- Evidently, the three operators 
X, y, z also satisfy the angular-momentum algebra, but 
with an effective Planck constant ?ieff = 1/ J. Taking the 
^off limit with fixed 77 and a, the classical dynamics 
associated with F can be obtained, with variables x, y. 
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FIG. 1: The eigen-phase spectrum (denoted e) of the Floquet 
operator F in Eq. (3). J = 20 in (a), 30 in (b), and 30.5 
in (c). a/TicS = 1 in all panels. Beautiful butterfly patterns 
are clearly recognized. Analogous result for J = 100 can be 
found in Ref. 19]. 
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FIG. 2; The dissolving of the butterfly pattern of the eigen- 
phase spectrum of the Floquet operator F in Eq. (3). In all 
the panels J = 30. From top to bottom: a/TicS = 3, 6, and 
10, respectively. This process resembles that observed in the 
kicked Harper model Q. 



and z restricted on a unit sphere. Because rj = AgoN , 
this classical limit with fixed t] requires N — > +oo and 
(7o — > 0. This condition is apparently equivalent to that 
in a standard mean-field limit of the driven BEC. 

In addition to the system defined by F, we also con- 
sider one of its interesting variants: 

Fxy = e''''^^/'^'^^e~*"'^="e~*''"'^/*^^'^''e~*"'^'', (6) 

which is different from F in the last factor, i.e., e"*"''^ 
in F is replaced by e"*"'^" . As seen below, such a variant 
may induce considerable changes in the spectral proper- 
ties. 



III. DETAILED ASPECTS OF THE 
BUTTERFLY SPECTRUM 

A. Multifractal properties 

In the I to) representation, the matrix elements of the 
operator F can be evaluated in a straightforward man- 
ner. Diagonalizing F numerically then yields its spec- 
trum. Figure [1] shows the typical eigen-phase spectrum 
of F YS hrf = "qhee = "q/J = 8go, for J = 20, 30, 30.5 and 
ct/^oS = ctJ = 1- Because the spectrum of F is invari- 
ant if hri + 47r (see the proof in Appendix A), we 
set hr^ € [0,47r). Though in Fig. [T]the involved Hilbert 
space is rather small, spectacular butterfly patterns are 



already obtained. Their reflection symmetry with respect 
to h,i = 2-K is also clearly seen, a fact proved in Appendix 
B. The found butterfly patterns in Fig. 1 resemble the fa- 
mous Hofstadter's butterfly, but also present remarkable 
differences in several aspects. First, if we take a vertical 
cut of the butterfly patterns in Fig. [1] the spectrum is 
not found to present any large gaps. Second, the but- 
terfly patterns shown in each panel of Fig. [1] possess a 
double-butterfly structure, with each butterfly covering a 
27r range of firi- This double-butterfly structure is some- 
what analogous to the spectrum of a Harper-like effective 
Hamiltonian considered in Ref. [ll[. More interestingly, 
though Fig. mb)-(c) has more levels than Fig. [Tfa), the 
overall outline of the double-butterfly structure is seen to 
be insensitive to J for fixed a/hcs = aJ. Indeed, in Fig. 
1(c) of Ref. [l^ we also presented the spectrum for a 
much larger J value, i.e., J = 100, and again similar out- 
line of the butterfly spectrum is obtained. Qualitatively 
this is because when a/hcs ~ aJ is fixed, the phase 
range of the second and fourth factors of F is also fixed. 
By contrast, for a fixed value of J but for other not too 
large values of a, the qualitative features of the butterfly 
spectrum remain, but at different scales. For very large 
values of a (e.g., a/?icff > 10), the butterfly pattern for a 
fixed value of J will gradually dissolve, as seen in Fig. [2l 
This dissolving process of a butterfly spectrum is similar 
to that seen in the kickcd-Harper model 

Some detailed features of the spectrum are also note- 
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FIG. 3: Cumulative Floquet state density (a) and the Floquet 
state density distribution (b)-(d) at different scales, for = 
(75 - l)7r/2, a/TicB = 1, and J = 2999. Panel (e) shows 
the generalized fractal dimension Dq. Crosses and circles are 
for odd-parity and even-parity states. Triangles represent the 
result for a standard kicked-top model. 



worthy. For example, it is observed that the spectrum 
collapses to one point for hri — 27r, if and only if J is 
an integer. This can be explained as follows. If J is an 
integer and if hjj = 27r, then in the \m) representation, 



„ — ITT 



(7) 



So in this case e^"''^^/'^-^^ is equivalent to a rotation of tt 
around the z axis, and hence the first three factors of F 
exactly cancel its last factor. This cancellation will not 
occur if J is a half integer, i.e., if N is odd. Later we will 
return to this intriguing difference between odd-A'^ and 
even-iV cases. 

We have also examined the statistical behavior of the 
found butterfly spectrum. To have good statistics we 
consider a much larger J = 2999. Figure [31Ja) presents 
the cumulative level density iV(e) for a representative 
value of hri- It is seen that N{e) is highly irregular, but 
does not show any clear flat steps. This is consistent 
with our early observation that no large gap exists in the 
spectrum. Figure [3fb)-(d) shows the associated level dis- 
tribution P(e) at three different scales. Evidently, P(e) 
has a fascinating self-similar property. This motivates 
us to quantitatively characterize the spectrum via the 
generalized fractal dimension Dq, with the results shown 
in Fig. [3{e). As expected from the N{e) result in Fig. 
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FIG. 4: (a) Eigen-phase spectrum of the Floquet operator 
Fxy (Eq. for J = 30. (b) The corresponding fractal 

dimension Dq computed for J — 2999 and fi^ = (\/5— l)'''"/2. 
a /heft = 1 in both panels. 



[3l[a), Do = 1. However Dg for g ^ clearly shows that 
the spectrum has multifractal properties. For the sake 
of comparison. Fig. ^e) also shows the Dq result for a 
standard kicked-top model with the same values of r/ and 
a (i.e., considering an operator comprising only the first 
two factors of F) . The Dq behavior in the kicked-top case 
is as trivial as that of a random sequence: it remains close 
to unity and slightly decreases with increasing q due to 
finite-size effects. Based on these results, we conjecture 
and invite a formal mathematical proof that the butter- 
fly patterns found here contain true fractals in the limit 
of J ^ -l-oo. 

We have also studied the F^y model deflned above, 
obtaining a similar multifractal butterfly spectrum, as 
shown in Fig. IH^a). Interestingly, despite that the out- 
line of the butterfly spectrum of F^y is much similar to 
that for F, careful investigations reveal considerable dif- 
ferences between the butterfly spectrum of F^y and that 
for F. For example, the fractal dimension Dq shown in 
Fig. [H^b) for Fxy is similar to, but slightly larger than, 
that of F for q > 0. It is found that this is because the 
gaps in the butterfly spectrum of F^y is more densely 
filled than that of F (compare Fig. g] (a) and Fig. W[h)). 
In the next subsection we will point out an even more 
fundamental difference between these two systems. 



B. Level crossings 

In this subsection we study the level crossings in the 
butterfly spectrum as the parameter varies. Note first 
that due to the above-mentioned symmetry (rn|_F|n) = 
{—m\F\ — n), the eigenstates of F can be classified into 
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J eigenstates of odd-parity and J + 1 states of even par- 
ity. As such, we should investigate the crossings between 
different-parity states and between same-parity states. 
In either case, computationally it is found that the mini- 
mal distance in hrj between two level crossings decreases 
sharply with J. So even for a rather small J ~ 10 it is 
already numerically demanding to identify all the level 
crossings. 

As an example Fig. [5lja)-(b) presents the typical level 
crossing behavior for J = 10. The Floquet states are 
seen to cross each other frequently, between different- 
parity states and between same-parity states. Both types 
of level-crossings turn out to be of vast interest. For 
the first type, at a crossing point an arbitrary superposi- 
tion of two crossing states of different parities remains an 
eigenstate but generally breaks the parity symmetry. So 
if such a superposition state is used as the initial state, 
the ensuing dynamics will maintain a nonzero popula- 
tion difference between the two modes forever [15[ . This 
makes it clear that the first type of level crossings give rise 
to the seminal CDT phenomenon [l^ that has attracted 
broad experimental and theoretical interests. It should 
be pointed out that in some regimes of to the naked 
eye two curves of opposite parities in Fig. [Sfa) and Fig. 
[5]Jb) are almost on top of each other. As a result many 
CDT points are found in these regimes. Note also that 
the CDT-induced population trapping is fundamentally 
different from the well-known self-trapping effect on the 
mean-field level. Indeed, the CDT effect here depends 
on 77 and J, whereas mean-field self-trapping is transient 
and independent of J. 

Now turning to the second type of level crossings, they 
come as a surprise because avoided crossings between 
same-parity states, rather than true level crossings, are 
generally anticipated for classically non-integrable sys- 
tems (see Figs. HHT]). The second type of crossings 
therefore suggest the uniqueness (e.g., some effective lo- 
cal "symmetry") of F whose matrix elements in the |m) 
representation are quasi-periodic. Recalling the above- 
mentioned extreme example where all levels cross at 
= 2tt for integer J, we expect that special arithmetic 
properties of play a key role in both types of level 
crossings. 

Careful checks are made to ensure that the same-parity 
level crossings observed here are not avoided crossings 
with a very small gap. For example, we examined the 
crossing behavior for small J, where analytical studies 
become possible. In particular, for J = 2, using Wigner's 
rotational matrices to express the 2nd and 4th factors of 
_F, we can analytically prove that there must be true level 
crossings between two odd-parity states, at = 2tt/3 
and = IOtt/S, regardless of the value of a. This is fully 
consistent with our numerical finding. Details for this 
case are presented in Appendix C. This further confirms 
that the number theory properties of are responsible 
for the same-parity level crossings. As another check, 
we also studied the level crossing behavior in the but- 
terfly spectrum of F^y. Therein we only obtain avoided 
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FIG. 5: (color online) (a) Level crossings between 10 even- 
parity states (solid) and 9 odd-parity states (dotted) for J = 
10 and a/hcd — 1.0. (b) A magnification of one part of (a), 
(c) The number of level crossings versus J, for hrj £ [0, 4n) 
and a/TicS ~ 1.0. The cross (square) symbols are for crossings 
between different-parity (same-parity) states and the fitting 
suggests a power law scaling ° (J'^'^). 



crossings between all the eigen-phases (see also Fig. 9 
in Appendix C). For J < 12 and for the same parame- 
ters as in Fig. OJc), the typical gaps of the avoided level 
crossings in the spectrum of F^y are found to be > 10^^, 
many orders of magnitude larger than the accuracy of 
our eigen-phase calculations (10~^^). This "control case" 
hence indirectly supports our observation of same-parity 
level crossings for F. 

By obtaining all the level crossings in the butterfly 
spectrum of F with high accuracy for J < 12, we obtain 
in Fig. OJc) that the number of CDT points contained in 
the butterfly patterns scales as J'^ " and the number of 
same-parity crossings scales as J^ "^. In either case, the 
number of crossings divided by the total number of levels 
(~ J) or divided by the total number of level pairs (~ J^) 
diverges as J — > +00. In particular, we assert that as N 
goes to infinity, on average each pair of Floquet states in 
a butterfly pattern see infinite CDT points. 



C. Relevance of the classical limit 

An interesting question is what implications the classi- 
cal dynamics may have on the fractal spectrum observed 
in the F system. To this end we investigate the classical 
limit of F defined above. Numerical studies indicate that 



6 



?S 



?N0 




?\ 



?S0 




FIG. 6; Poincare surfaces of section (with Jx > 0) of the 
classical or mean-field limit of F in Eq. (3), with a = 0.05 
(same as in Fig. 1(a)), 77 = 5 in (a), 30 in (b), 75 in (c) and 
100 in (d). 



classically the system can be governed by both regular 
and chaotic motions, depending on the two system pa- 
rameters 77 and a. In general, for a fixed a, as 77 increases 
the classical dynamics undergoes a transition from being 
regular to being chaotic. As an example Fig. [S] presents 
the phase space structure of the classical limit of F with 
a = 0.05 and an increasing rj. 

Such a classical regular-to-chaos transition lacks a 
quantum counterpart in the butterfly spectrum shown in 
Fig. [l]and Fig. [2l whose characteristics can be much sim- 
ilar for radically different values of rj. Indeed, the quan- 
tum Floquet spectrum is periodic in 77 with a period 4 Jtt. 
(See Appendix A.) Therefore, upon quantization the reg- 
ular or chaotic nature of the classical dynamics may not 
necessarily be reflected in the spectrum and hence can 
be irrelevant to the quantum dynamics. 

Interestingly, the regular-to-chaos transition in the 
classical dynamics is not even connected to the dissolving 
of the butterfly spectrum (see Fig. [2|). As an example, 
let us examine such a dissolving process for J = 30 from 
a/h^s = 1 [Fig. [2b)] to a/h^s = 3, 6 and 10 (Fig. 

corresponding to a = 0.1/3,0.1,0.2, and 1/3 respec- 
tively. Consider the spectrum along the vertical cuts of 
these figures at Ti^/tt = 2 x 10"'^, i.e., 77 = O.OGtt. It is 
found that in spite of the drastic changes of the spectrum 
in this dissolving process, the dynamics in the classical 
limit remains regular all the way from a = 0.1/3 to 1/3 
[see Fig. [7Ua)-(b)]. On the other hand, due to the 47r 
periodicity of the quantum spectrum in we expect 
the same dissolving process if we change 77 from O.OGtt 
to O.OGtt -t- 1207r. However, the corresponding classical 
dynamics is now completely chaotic [see Fig. [7Uc)-(d)]. 
Hence the dissolving of the butterfly spectrum is unre- 
lated to the classical dynamics. This is in contrast to 
the kicked-Harper model, where the dissolving of Hosf- 



FIG. 7: Same as in Fig. [6]but a — 0.1/3 in (a) and a = 1/3 in 
(b) with rj = O.OGtt; a = 0.1/3 in (c) and a = 1/3 in (d) with 
Tj — O.OGtt -1- 120tt. Quantum mechanically, for J — 30, cases 
(a) and (c), or cases (b) and (d), share the same spectrum 
due to the spectral symmetry under Tirj Tin + 4tt. 



tadter's butterfly was somewhat connected to the classi- 
cal regular-to-chaos transition [5] . 

The little relevance of the classical limit in under- 
standing the butterfly spectrum also makes our model 
markedly different from the conventional kicked-top 
model as a paradigm for quantum chaos. We believe 
that such a lack of classical-quantum correspondence is 
a result of a built-in feature of our SU(2) model. That 
is, our strategy for generating the fractal spectrum is 
based on a partial cancelation of quasi-random dynam- 
ical phases (see Eq. [S]), and this partial cancelation of 
quantum phases is a pure quantum feature with no clas- 
sical analog. To double check this we also studied the 
classical dynamics of the F^y model and reached the same 
conclusion. 



D. Generating a fractal-spectrum family 

In the previous subsection we have studied the but- 
terfly spectrum associated with F . Here we point out 
that the butterfly spectrum of F is just one member of 
a whole butterfly-spectrum family. Let us restrict our- 
selves to cases of integer J and consider the following 
function 

9{t) = 5I[5o^(< - 271T - r) + ~g^5{t - 2nr)] 

n 

instead of that used in generating F as in Eq. Here 
go (§ ■ ^ — 170) and ji and v are two integers sharing no 
common factors. We can then obtain an extended class 
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of Floquet operators 

We also note that i''^''/'') can be obtained if we set 
45ot/'C — 2771^//^ in Eq. H]). Obviously, F corresponds 
to the special case of = /i = 1. 

In the double-kicked-rotor model considered in Ref. 
01, one can construct analogous operators by employing 
high-order quantum resonances, where v / ^, indicates the 
resonance order. For each choice of v / a certain type 
of fractal spectrum can be generated. An example for 
V / fjL = 1/2, the so-called anti-resonance condition, can 
be found in Ref. 0| (see Fig. 3 therein). 

Interestingly, the spectrum of F*^"/^' defined in Eq. 
([8|) also forms a fractal-spectrum family, with the out- 
line of each butterfly roughly similar to its relative in the 
double-kicked-rotor model considered in Ref. [7j . We also 
find that for integer J, the peculiar spectral characteris- 
tics of F discussed above can be maintained in the spec- 
trum of F^'^/'^K For example, we have investigated thor- 
oughly the case of i^/fi = 1/2 and obtained that (l) Dq 
curve is qualitatively the same as that for F; i.e. Dq = 1 
and Dq decreases as q increases, (2) level crossings be- 
tween eigenstates of the same or different parities can 
occur; and (3) the dynamics in the classical limit is also 
irrelevant in understanding the quantum spectrum. Sim- 
ilarly, the spectral properties of F^y are also found to be 

similar to its high-order extension Fxy^^^ = e^"^^ / Fxy ■ 
The possibility of constructing such a butterfly- 
spectrum family provides further support that our strat- 
egy for generating a fractal Floquet spectrum in driven 
quantum systems is quite general. 



IV. DISCUSSION 

Experimental confirmation of a butterfly Floquet spec- 
trum in driven systems is challenging. In the case of the 
kicked-Harper model, there have been a few experimen- 
tal proposals but so far the kicked-Harper model has not 
been experimentally realized. The double-kicked-rotor 
model proposed in Ref. Q opens up a new opportunity. 
However, in atom-optics realizations of the kicked-rotor 
model, a dilute cold gas in a kicking one-dimensional 
optical-lattice potential is required, the quasi-momentum 
spread of the initial state should be sufficiently narrow, 
and the interaction between the atoms should be negligi- 
ble. By contrast, in the present study, we rely on an in- 
teracting cold gas distributed on two modes, and there is 
no quasi-momentum issue. Moreover, because the effec- 
tive Planck constant is simply given by 850 j by tuning 
the atom-atom interaction constant alone we may scan 
the butterfly spectrum already. These advantages make 
the systems proposed in this study one more step closer 
to possible experiments of a butterfly Floquet spectrum. 



There are also other motivations for further theoret- 
ical and experimental studies of our model. First, be- 
cause the found butterfly spectrum collapses at = 2tt 
(or go — i'/4) for integer J, one may experimentally de- 
termine if N is even or odd by scanning the dynamics 
in the neighborhood of go = Note that detecting 

the even-odd properties of TV is impossible in the mean- 
field dynamics of a BEG, and such a topic is already 
under investigation in Ref. [25j using a different mecha- 
nism. Similarly, one may study the CDT points to reveal 
non-mean-field effects. Second, it is now of great inter- 
est, both experimentally and computationally, to revisit 
early results of how a multifractal spectrum can be mani- 
fested in time-dependent properties [2^ . Third, noticing 
that recently dissipative two-mode BEC's have attracted 
considerable attention [23|, it seems interesting to study 
dissipation effects on a butterfly spectrum. 

Experiments of our model system with a small J can 
be also interesting. For example, for the J = 2 case we 
have found that the F operator is an identity matrix in 
the odd-parity subspace at Ti,^ — 27r/3 and 107r/3. (See 
Appendix C). Such cases may be used to extract the value 
of go to a good precision. They may be also useful to 
explore the implication of level crossings between same- 
parity states. 

We next discuss how to experimentally confirm a but- 
terfly spectrum. Certainly, a Floquet spectrum cannot be 
directly measured. However, in principle there is a stan- 
dard procedure to invert the Floquet spectrum from the 
dynamics. This involves quantum state reconstruction of 
the evolving two-mode BEC. That is, after each period of 
the driving field, enough measurements on different ob- 
servables are to be made to reconstruct the state of the 
evolving two- mode BEC. This is possible because quan- 
tum state reconstruction of a BEC is currently a very 
active and fruitful area. This state reconstruction is not 
expected to be prohibitively demanding if N is relatively 
small. With the quantum states at different times recon- 
structed, then by making a Fast-Fourier- Transformation 
(FFT) of the time-evolving state, the Floquet spectrum 
may be obtained. 

We point out that there is no need to follow the dynam- 
ics for very long in order to resolve a butterfly spectrum. 
Let us assume that the wavefunctions are already recon- 
structed by sufficient measurements at integer multiples 
of the driving period r. To examine how many periods 
are needed for obtaining the butterfly spectrum from ex- 
periments, we consider an example where the initial state 
is chosen as |m = 10) for J ~ 20. Using the same param- 
eters as in Fig. [ija), and using the time-evolving states 
after n = 32, 64, 128, 256, 512 and 1024 kick pairs, the 
spectrum obtained via FFT are shown in Fig. [Sj It is 
seen that a few hundred kicks can be good enough to re- 
solve the shape of a butterfly spectrum relatively well. If, 
as discussed in the following, the time scale of the kicking 
period is chosen to be 10~^s, then this means that the re- 
quired duration to follow the dynamics in an experiment 
is around 0.01s. 
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FIG. 8: Eigen-phase spectrum of the Floquet operator F for J = 20 and a/TicH — 1, retrieved via an FFT of the sequence 
{((/>(O)|F"|0(O)),n = I,-- - ,N} with the sequence length n = 32(a), 64(&), 128(c), 256(d), 512(e) and 1024(/). The initial 
condition is given by |<?!'(0)) — |10) as an example. 



Finally, let us comment on the parameter ranges we 
have chosen. First of all, to ensure that a wide regime 
of the found butterfly spectrum can be visited in a real 
system, a tunable go is required, and its characteristic 
value should be ~ 1. Let gc be the self-interaction con- 
stant of a static two-mode BEC. A reasonable range of 
9c/ f if is the tunneling rate) is from 10"'^ to 10~^ for a 
two- mode BEC in a double- well potential [13]. With the 
first realization of F in mind and for a J = 1 and J = 10, 
we have gc/ f = '^^Tg^/a, and that 2Tgc ranges from 10^^ 
to 10"'^. This indicates that go/irgc) should be around 
10'^ — 10^. If the value of gc (in SI unit of frequency) is 
about 50 s'-^ (a value considered in [2Q|), we have that 
r is in the range of 10~^s to 10~^s. Considering other 
realizations of a two-mode BEC might lead to different 
characteristic values of r and go/irgc)- 



V. CONCLUSION 

We have presented a strategy for generating a fractal 
butterfly Floquet spectrum in a class of driven SU(2) 
systems. The essence of this strategy is to partially 
cancel the quasi-random dynamical phases in the time 
evolution and then induce intriguing spectral properties. 
The success of such a strategy in both this work and 
our early work treating a double-kicked-rotor model 
indicates its wide applicability. As such, butterfly Flo- 
quet spectrum is expected to occur in many driven quan- 



tum systems that can go well beyond the context of two- 
dimensional electron systems or the Harper model con- 
text. A butterfly spectrum thus obtained may also differ 
significantly from Hofstadter's butterfly. 

Detailed aspects of the butterfly spectrum in driven 
SU(2) systems are also examined. The level crossing fea- 
tures in the butterfly spectrum, especially its connection 
with coherent destruction of tunneling, and the surpris- 
ing observation of level crossings between same-parity 
states, are emphasized. The sensitivity of the butter- 
fly spectrum on the even-odd properties of the number 
of particles is also highlighted. Many further questions 
can be asked and we believe that our driven SU(2) model 
will be relevant to a number of related research topics, 
including fractal statistics, quantum chaos, quantum con- 
trol, and the validity of mean-field descriptions of Bose- 
Einstein condensates. 

The conventional system for understanding a butterfly 
Floquet spectrum is the kicked-Harper model. Therein 
the quantization rule varies with the boundary condition 
adopted ^2§\ and a compact toroidal phase space arises 
only if the Planck constant assumes special values ^29j |. 
A general treatment of the kicked-Harper model leads to 
a band structure that often complicates the issue. By 
contrast, the phase space structure in our driven SU(2) 
model is necessarily on a sphere j2^ . with no arbitrari- 
ness in quantization and no band structure in the spec- 
trum. For these reasons we hope that the butterfly Flo- 
quet spectrum discovered here can stimulate future stud- 
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ies on general implications of a fractal Floquet spectrum. 
Combining this work with our early study we have 
that three paradigms of quantum and classical chaos, i.e., 
the kicked-rotor model, the kicked-Harper model, and the 
kicked-top model, are linked all together because any one 
of the three contexts can be used to generate quantum 
critical systems with a fractal Floquet spectrum. 
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APPENDIX A: 4tt PERIODICITY OF F IN ?i„ 



Denoting by U{hjj) = e'''^^/(^'^\ the operator F can be 



written as F{Tiri) 
shown that 



C/(?i^)e-'"-^-C/t(;i^)e 



C/(?i^ + 47r) = ±U{hr,), 



It can be 



(Al) 



hence + 47r) = To show the above equality 

let us consider the representation using the eigenstates 
of Jz, i.e. Jzim) = nih\m)^ in which U is diagonalized 
and {m\U{%^ + 47r)|TO) = ^lin+'iJ^W /(^J) ^ resulting in 
U{fir) + 47r) = U [fin) if J is an integer, and U(fir) + 47r) — 
— U{h,^) if J is a half- integer. In the latter case as {hrf + 
47r) = —U^hjf), which contributes a second minus sign, 
we again have F{hji + in) = F{hjf). 



APPENDIX B: SPECTRUM SYMMETRY OF F 



UNDER h„ 



-hr, 



Consider first the F operator and its eigenfunction Itjj) 
with 



(Bl) 



Multiplying e ^KJ^/^g ""^J^^ to both sides of the above 
eigenfunction equation and defining that = 

g-in^j^/2g-iQja: 1^^^ ^j-^gj^ have 



'W) 



(B2) 



Hence, comparing Eqs. (jBip and (jB2[) one sees that the 
spectrum of F{hri) is identical with that of F{—fi,-i). 

Combining this result with the Att periodicity of F in 
firi, one obtains that the butterfly spectrum of F has a 
reflection symmetry with respect to fin — 27r. 



0.4-1 
^ 0.0^ 
-0.4 




0.4-1 
Uj 0.0- 
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FIG. 9: (color online) The eigen-phase spectrum of the Flo- 
quet operator F (a) and F^^y (b) for J = 2 and a/hcft — 1. 
In (a) the two eigen-phases of odd parity (dotted) cross at 
hr, — 2tt/3, 2tt, and IOtt/S, while the three eigen-phases of 
even parity (solid) only cross at hn = 2-k. As a contrast, the 
spectrum of F^y displays avoided crossings instead. 



APPENDIX C: LEVEL CROSSINGS BETWEEN 
TWO ODD-PARITY STATES OF F FOR J = 2 

The arithmetic properties of hrj can cause surprising 
and frequent level crossings between same-parity eigen- 
states of F. To understand this counter-intuitive result 
we consider a simple case with J = 2. Numerically we 
observe level crossings between two odd-parity states at 
h,, — 27r/3 and 11^ — 107r/3 besides that at — 27r. 
(see Fig. 9a for a/Tics — I)- Here we give a mathemat- 
ical proof that this is indeed the case, thus supporting 
the numerically observed level crossings in general. Note 
that our approach can be extended to other cases as well. 

We first note that each of the four factors of F pre- 
serves the parity. So if we consider the odd-parity sub- 
space only, then each of the four factors can be reduced 
to a 2 X 2 matrix because the odd-parity subspace is 
two-dimensional for J — 2. Specifically, we consider the 
following basis states 

|2) ^ (|2)-|-2))/a/2 
11) = (ID -|-1))/V2 



In this representation, the factor e^^^'-^^^^ and e ^f^vhf^ 
for — 27r/3 are given by 



U : 



gi47r/3 

e*''/3 



and respectively. 

In order to get the matrix expression of F, we need to 
get the analytical expression for the factor e~'""'^ . Fortu- 
nately, this can also be done by using Wigner's rotational 
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matrices Denoting -D^/,„ as Wigner's expression for 
the rotational matrix, we have 

Dl,„,^{rn'\e-^'^-'y\m). (CI) 

What we need is the matrix (m'|e~*"'^^ |to), which is re- 
lated to f^j/^ by 

■^^|m)=z"'-'»i?,V„. (C2) 

Using the explicit expressions of the rotational matrix 
^m'm given in [3Q| , we obtain 

(2|e^*"^-|2) = cos4(a/2) 

(2|e"'"^-| -2) = sin4(a/2) 

(l|e-*"^-|l) = [l + cos(a)][2cos(a)-l]/2 

(l|e-*"^-| - 1) = -[l + 2cos(a)][l -cos(a)]/2 

(2|e-*"'^-|l> = -isin(a)[l + cos(a)]/2 

(2|e-'"^-| - 1) = -isin(a)(cos(a) - l)/2 

(l|e-*"^-|2) = -isin(a)[cos(a) + l)]/2 

(lle"*"-^-! - 2) = -isin(a)[cos(a- l)]/2. 

In the odd-parity subspace, we then find 

(216-'"-^== |2) = cos(a) 

(216-'"-^- 1 i) = -2sin(a) 

(i|e-'"'^-|2) = -zsin(a) 

(ile-'^-^-li) = cos(a). 



Finally, the i^-matrix in the odd-parity subspace is given 
by 



Because F turns out to be a unit matrix, this directly 
demonstrates that the two eigenphases of F are both zero 
(because — 1). Two odd-parity eigenstates hence 
cross each other at — 27r/3. The symmetry of the 
spectrum then indicates another crossing at = IOtt/S. 
In addition, this proof shows that this crossing is inde- 
pendent of a, which is also confirmed by our numerical 
calculations. Our proof can also be extended to other 
cases with a rather small J, because thanks to Wigner, 
the rotational matrix elements can be analytically ob- 
tained. 
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